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Abstract 

Some general aspects of nonlinear transport phenomena are discussed on the 
basis of two kinds of formulations obtained by extending Kubo's perturba- 
tional scheme of the density matrix and Zubarev's non-equilibrium statistical 
operator formulation. Both formulations are extended up to infinite order of 
an external force in compact forms and their relationship is clarified through 
a direct transformation. In order to make it possible to apply these formu- 
lations straight-forwardly to thermal disturbance, a mechanical formulation 
of it is given (in a more convenient form than Luttinger's formulation) by 
introducing the concept of a thermal field which corresponds to the tem- 
perature gradient and by defining its conjugate heat operator Ah = Ylij ^j'^'j 
for a local internal energy hj of the thermal particle j. This yields a trans- 
parent derivation of the thermal conductivity k of the Kubo form and the 
entropy production {dS / dt)-^-^- = kE^/T. Mathematical aspects on the non- 
equilibrium density-matrix will also be discussed. In Paper I (Physica A 
390(2011)1904), the symmetry-separated von Neumann equation with re- 
laxation terms extracting generated heat outside the system was introduced 
to describe the steady state of the system. In this formulation of the steady 
state, the internal energy (Ho)f is time-independent but the field energy 
(T^i)^ (= — {A)^- F) decreases as time t increases. To overcome this problem, 
such a statistical mechanical formulation is proposed here as includes energy 
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supply to the system from outside by extending the symmetry-separated von 
Neumann equation given in Paper I. This yields a general theory based on 
the density-matrix formulation on a steady state with energy supply inside 
and heat extraction outside and consequently with both (Ho)^ and con- 
stant. Furthermore, this steady state gives a positive entropy production. 

The present general formulation of the current yields a compact expres- 
sion of the time derivative of entropy production, which yields the plausible 
justification of the principle of minimum entropy production in the steady 
state even for nonlinear responses. 

Keywords: irreversibility, entropy production, principle of minimum 
entropy production, transport phenomena, electric conduction, thermal 
conduction, linear response, energy supply, steady state, Kubo formula, 
symmetry-separated von Neumann equation, Zubarev's non-equilibrium 
statistical operator. 



1. Introduction 

In the first paper of this series [l[ , the present author has derived the irre- 
versible entropy production in transport phenomena [1-16] starting from the 
von Neumann equation, as Kubo et al. [2-4] formulated the linear response 
scheme. A new aspect of the present author's theory [ij is that the entropy 
operator S is defined using the equilibrium density matrix peq as 

S = -kB logpeq = {'Ho - Fq)/T (l) 

with Fq = — /cBTlogTrexp(— /37/o) and that the entropy production is ob- 
tained thereby from the symmetric part Ps{t) of the density matrix[l|: 

Here, Hq denotes the Hamiltonian of internal energy, F a static field, ap the 
nonlinear transport coefficient and T the temperature. Thus, the dissipation 
(or Joule heat in electric conduction) apF'^ is obtained [ij from the irreversible 
entropy production as apF"^ = T{dS /dt)[„ even in the nonlinear regime 
as well as in the linear scheme. This is a big contrast to the traditional 
phenomenological formulation of dissipation using the complex admittance. 
On the other hand, the current Jp is given hj Jp = Trjtpa(^) = crpF using 
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the antisymmetric part p^it) and in particular the Kubo formula is given by 
the first-order term pi{t) of the antisymmetric part Pa{t) as 

POO pP 

J = Tijpi{t) = cJoF; ao= dt dXe~'' {jjit + ihX))^ (3) 

Jo Jo 

with the adiabatic factor e~^*. Here, the density matrix p{t) is separated 
more explicitly into the two parts ps(t) and Pa(^) defined by 

Ps{t) = po + P2{t) + Pi{t) + ■ ■ ■ , (4a) 

and 

Pa(t)=Pl(t)+P3(t)+P5(t) + ---, (4b) 

as the power-series expansion with respect to the field F{t) defined in the 
Hamiltonian 

H{t)=Ho + Hi{t)] Hi{t) = -A-F{t). (5) 

One of the purposes of the present paper is to apply the general theory 
[T] to the thermal conduction by constructing a mechanical formulation of 
the thermal conduction. For this purpose, the concept of a thermal field 
and its conjugate heat operator Ah will be introduced in Section 5. This 
formulation yields a new direct derivation of the heat current operator jf ^ in 
the Kubo formula ([3]). 

Before presenting the mechanical formulation of thermal conductivity, 
it will be instructive to formulate nonlinear responses in two ways, namely 
using the Kubo-type scheme and by extending the Zubarev-type scheme up 
to infinite order of the field F{t). The relationship between these schemes 
will be studied. The unboundness of the relevant non-equilibrium operators 
obtained in the present paper will be also discussed. 

Finally the idea of the symmetry separation in the von Neumann equa- 
tion is shown to be useful in formulating the energy supply to the relevant 
system, without changing the previous formulation pj of the internal energy 
('Ho)t, current Jp and irreversible entropy production {dS/dt)iT:r- For this 
purpose, we introduce a new concept of the stationary density-matrix whose 
time derivative can be allowed to be not only zero but also a non-vanishing 
constant operator. 
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2. Kubo-type formulation of the non-equilibrium density matrix 



In the previous paper the author studied the irreversibility and en- 
tropy production in transport phenomena [2-16], using the idea of the sym- 
metry separation of the density matrix in the von Neumann equation. 

It will be instructive to discuss explicitly the nonlinear solution of the 
von Neumann equation 

^o{t) = hH{t),p{t)]- U{t) = -Ho + n^it) and Ui{t) = -A ■ F{t) (6) 
at iTi 

for a field F{t). Following KuboQ, we may put p(t) = Po + ^p{t) with 

po = e-^^VZo(/3); Zo(/3) = Tre'^^". (7) 
The solution Ap{t) satisfying the equation 

|Ap(t) = ^[n{t),Ap{t)] + i-[Hi(t),po] (8) 
is given formally by 

= 4 fmt'm,{t'),po]i^-\t,t')dt', (9) 

where U(t,t') is the ordered exponential defined by [17-21] 

U{t,t') = exp+ (^^^^ n{s)ds^ 

I ft / X \ ^ /•* r^i 



^ (^) / '^^^ It' ' " /' '^^"'^(^1)^(^2) ■ • • Hitn) + 



(10) 



It will be more convenient for expressing the current Jp = Trjp(t) = 
TrAAp(t) in the Kubo-type (generalized) form 22|, |23| to rewrite Ap(t) as 

(Ap(t))K-type = ^ dt'U{t,t') (^j^Xpoji-inX)^ U~\t,t') ■ Fit'), (11) 
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where 

Q{t) = e-^^''Qe^^° and j{-ihX) = e^^«je-^^0; j = A= ^[A^o] 

(12) 

for any operator Q excluding the density matrix. Thus, the nonhnear current 
Jp is given in the form 

Jp = Trj (Ap(t))K.t,pe = fdt' f d\ {3{-in\)U-\t,t')jU{t,t'))^ ■ F{t'). 

J to Jo 

(13) 

If the force F{t) contains only the adiabatic factor e'^*(e — )■ +0), namely 
it is essentially time-independent as in Paper I, then the above expression 
ffTTj) is simplified as 



(Ap(t)) 

K-typc 



J^^ dt' (^j^^ dXpoji-thX)^ it' -t;F)- Fe'''. (14) 
Here, Q{t; F) except the density matrix p{t) is defined by 

Q(t; F) = exp (" J"^) Q exp J-H^ ; Q(t; 0) = Q{t) (15) 

for "H = "Ho + "Hi = "Ho — A ■ F and consequently we obtain Jp = apF, 
where 

/•oo pB 

ap= dte~'' dX{j{-ihX)j{t-F))^ ■ 3{t-F) = e-^je^. (16) 
Jo Jo 

Clearly this is reduced to the well-known Kubo formula ([3]), or more generally 

J = Re (ao(cj)Fe^'^*) ; (Jq{u) = dte^'^-'^' dX {jjit + ihX))^ (17) 

Jo Jo 

in the linear response scheme^], for the complex admittance cr{uj) calculated 
using the complex external field F{t) = Fe^'^^ in the general expression 
(fT3|) . As is well known, the first equation in Eq. (fT7|) is shown to agree with 
the expression directly obtained from Eq. (fT3|) for the hermitian Hamiltonian 
Hiit) = -A ■ Fcos{ut). 

This scheme is a statistical-mechanical perturbation theory, namely per- 
turbation with respect to the statistical ensemble, and it is not mechanical 16 
in contrast to the Zubarev-type formulation. 

The above formula f|T^ is useful in the plausible justification of the prin- 
ciple of minimum entropy production in the steady state, as will be discussed 
in Section 9. 
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3. Zubarev-type formulation 



There is anotlier type of formulation on nonlinear responses, namely the 
Zubarev-type non-equilibrium density-matrix scheme 0, S]- 

The starting point of the present general formulation is the following 
expression of the non-equilibrium density matrix 

p{t)=U{tM)pdA-\tM) = -^^exp{-f3U{t,to)'Hol(-\t,to)) 

expi-(3{'Ho + n{t,to))). (18) 



The non-equilibrium operator TZ{t,to) defined in Eq. ffTSl) satisfies the time- 
evolution equation 0] 

^7^(t,^o) = ^[n{t),nit,to)] + ^[niit),no], (i9) 

which is proved directly or using the quantum analysis [19-21]. It should be 
noted that the above equation ( |T9l) is purely mechanical, because it does not 
include the density matrix po (or pit)). The source term in Eq. (fT9|) is given 
in terms of the commutator ['Hi(t),'Ho] instead of the commutator [Tii, pq] in 
Eq.([8]). This corresponds to the Lippmann-Schwinger perturbation scheme 
in quantum scattering theory (j^l- 

Now the formal solution of 0191) is given by 

7^(t,to) = - [ dt'U{t,t')jU-\t,t') ■ F{t'). (20) 

Jto 

Then, the nonlinear current Jp is expressed by 

= Trjp(t) = Trj (Ap(t)),_,yp, , (21) 

where 

(Ap('^))z-type = Pit) -Po = -Po /'rfAe^^«7^(t,^o)e-^(««+^(*'*«)) 

Jo 

= -[ dAe-^(««+^(*•*"»7^(^,^o)e^^Vo, (22) 
Jo 
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or 



(Ao(t)) - ^ I du— 



d 

djji 

1 /-^ 



^o(/3) 



Jo 

(23) 



using the integration technique [231 . The expression (1221) is convenient for the 
perturbation expansion with respect to the field F{t), but it is asymmetric 
with respect to the A-integral. On the other hand, the triple integral rep- 
resentation (123|) is symmetric with respect to the field F{t) (namely to the 
/i- integral). Thus, the nonlinear current J p is expressed in the following two 
ways: 

Jp= f dt' rrfA(je-^(^°+^(*'*"%(t,OiW-^(t,t')e^''">oi^(t') 

J to Jo 

= f dt' d\{e^'^'Uit,t')jU-\t,t')e-^^^"+'^^'''"^^j)^Fit'), (24) 

Jto Jo 



or 



Jf = I dt' dX dfi 



to JO JO 

{e'^P~^^^'^^^+^'^i^M)U{t,t')jU-\t,t')^^^^ (25) 

as is easily seen from Eqs. (j22l) and (|23ll . respectively. The above expressions 
(IMl) and (!25ll yield immediately the transport coefficients, and particularly 
the new expression ( |25l) is interesting, because the quantum effect through 
the A-integral is taken into account symmetrically and gradually through the 
/i-integral. 

If the force F{t) depends on time t only through the adiabatic factor 
e'^*(e — )■ +0), then using the nonlinear current operator j{t]F) defined in 
Eq. flTBl) we obtain 

/»oo rl 

JF = crFF; aF= dte-'' dfi {j;j{t; F)) , (26) 
Jo Jo 
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in the limit to ~^ — oo, where 



= r dX (e-^^-m.) je-^^(-) jit; F)e^'^')^ (27a) 
Jo 



with 



and 



ni^) = no + fi hm 7^(^,to), (27b) 

to-^—oo 

n{t,to) = - c/t'e^* eT«je-T« . F. (27c) 

J to 

By the way, the Zubarev-type density-matrix is given in the form 

PZ-type = p(0) = exp i-PCHo + Uf)) /^o(/3), (28a) 

using 7^F = 7eF(0,-oo) = 7^i7(^, -oo), which shows W(t, 0)(?{o+^F)W-^(t, 0) 
Ho + T^-i? and consequently we have [Hq + 'R.f, H] = 0. (Note that this equa- 
tion can not determine TZp uniquely.) Here, from Eq. fl27cl) with Eq. f|T6l) . Tip 
is expressed by the integral 

POO 

nF = - e-''j{-t; F) ■ Fdt. (28b) 



Now we make here a linear approximation oiTZp with respect to the force F 
to obtain 



7^l = - / C^j{-t) ■ Fdt with j{t) = e^^^^je^^", 
Jo 



(29) 



if it exists in the limit to ~^ — oo. Then we obtain Zubarev's non-equilibrium 
statistical operator of the form 

Pzub=exp(-/3(?/o + 7^l))/Zo(/3) 

= exp (^-mo + /3 ^ e-^'ji-t) ■ Fdt^ /Zo{(3). (30) 

Note that [Ho + Tli,^] = 0(F'^). If we expand Eq. fl30|) again up to the linear 
order, then we obtain the Kubo formula (jSj). The above density matrices 
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(Ap) K-type and (Ap)z_type (or pz-type in Eq.(|28ap) are useful in studying a 
current operator, but we have to be careful in discussing the time change of 
physical operators such as Ho, Hi and entropy, in which the formula 

j^{Q), = ^Tr[Q,n]p{0) (31) 

should be used in the limit to —oo for any physical quantity Q, when the 
system is described by the original von Neumann equation ([6]) with 'H(t) = "H. 
The relation between the above Zubarev-type expression fl2^ and the Kubo- 
type expression (ITT]) will be discussed in the succeeding section. 

4. Relationship between the Kubo-type and Zubarev-type formu- 
lations 

As has been shown in Sections 2 and 3, the Kubo-type expression flTT]) is 
of convolution-type, and the Zubarev-type expression (!23|) looks very differ- 
ent from Eq. lITT]) . It will be instructive to prove the equivalence of the two 
expressions by a direct transformation. In order to realize clearly the proce- 
dure of the direct transformation, we explain the case of the time-independent 
non-equilibrium Hamiltonian Ti as follows. First note that 



(ApW)k..„. = (-P (^*«)) rn'^."' 



ih 



to 

t 



ew{lSn]]l dt' (A exp [-^n)) Po 



exp y—^^nj j Po- Po 

^ exp(-/3-Ho-Wt,to))-po = (Ap(t))z_t,pe. (32) 



ZoiP)^'^""- --v--^^ ru z_type 

where 5^ denotes the inner derivation[19-21] defined by = [H, Q] = 
HQ — QTi, and we have also used the relations (5^iPo = 5-hPo and 

exp {^-^^n^^ ^° " / (''''^ (^i^*^^)) ^[■^i''^o]+^o = n{t,to)+Ho. 

(33) 

This transformation is easily extended to the case of the general time- dependent 
field F{t). 
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5. Mechanical formulation of thermal disturbance to derive en- 
tropy production 

In order to apply the general formulations based on the von Neumann 
equation with the Hamiltonian form ([6]) presented in the preceding sections 
and to perform the first-principles derivation of entropy production even for 
thermal conduction, we have to formulate thermal disturbance in a mechan- 
ical way. 

The present theory [1] about entropy production is based on the von 
Neumann equation for density matrix p{t) and for the Hamiltonian 'H(t). 
The density matrix Ap(t) in Eq.(|H]) is statistical through pq, that is, it may 
contain the temperature T. Now we extend our interpretation of the von 
Neumann equation (jH]) so that we may treat thermal conductivity using the 
mechanical perturbation method in Kubo's paper jif. For this purpose, the 
Hamiltonian for thermal disturbance should be regarded to be an effective one 
renormalized in a local but enough large region to catch up the temperature 
gradient js], [oj, as in the Ginzburg-Landau- Wilson Hamiltonian describing 
critical phenomena. Then, we introduce a thermal field Et = V(3{r)/l3o as a 
field F{t) in Eq.([6]), where I3{r) = l/k-BT{r) and /3o denotes the characteristic 
/3. Such a renormalized particle as moves receiving a force from the field 
is called "a thermal particle". The next problem is to find an operator Ah 
conjugate to the thermal field Ei{t). The operator Ah should be chosen so 
that the time- derivative Ah may give the heat current Jh, namely Ah = Jh- 
The heat current Jh is first defined by 

•?H = ^{hjVj] = ^ J2^hjVj + Vjh,) (34) 

j 3 

in the simplest case of ideal renormalized particles, when the unperturbed 
Hamiltonian T-Lq is given in the form 

?/o = ^^^j; hj = £j+Uj with ~ 2 — ""j = ^(^i)- (35) 

3 

Here the last quantity Uj denotes the scattering by impurities, and m* denotes 
the effective mass defined by Eq. fl39|) . It should be also noted that [hj, hk] = 
or [efe,Mj] = for j ^ k. The time-derivative Ah is defined by 

An = ^[An,no], (36) 
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as usual. It is easily proved that the heat operator Ah of the form 



2 

satisfies the required relation 

Ah = ^[Ah,^o] = J2^h,v,} ^ jn- (38) 

j 

Here we have used the renormalized momentum pj and space coordinate Vj 
of the effective particle j defined in a local region of original n particles as 
follows: 

n 



Pj = Y.Pn^ ^J = ^Zl^i^' ^* = nm. (39) 

i=l 1=1 

We have also used the following commutation relations 

1 " 

[rj,Pj] = - ^[rji,pji] = ih and r-j = Vj, (40) 



n 

i=l 



which have been derived from Eg. (1391) . 

Thus, we have found a heat operator Ah conjugate to the thermal field 
in the form fl37|) . Once it is found, we may use it in a generalized case in which 
interactions among particles exist, that is, Uj is given by uj = u{rj, {r^}) 
and consequently we have [hj, hk] 7^ or [cj, m^] 7^ for such {A;} as the j-th 
and k-th particles interact. In this general case, the heat current is easily 
shown to be given in the form 



Now we expand p{t) in the von Neumann equation ([6]) as 

Pir(t) = Po + Pi{ty, Po = e-^^"/Zo{P) (42) 

with Zo(/3) = Tr exp(— /JHo), up to the first order of the thermal field E^. 
Then, the first-order density matrix pi{t) is expressed as 

p^(t) = Po f ds j dXET- Anis-t-ihX), (43) 

Jto Jo 
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where Ah denotes the current jjj of Ah, namely Jh = Ah = [Ah,?^o]/^^- 
Thus, we obtain the well-known formulajil, 0] 

J}i = {3n)t = Tr{(po + Pi(i)) Jh} = ^^^t, 

POO pfi 

K= ds d\e-''{j^j^{s + ih\))oET (44) 
Jo Jo 

by taking the limit to — )■ — oo for the initial to after the thermodynamic limit. 

The present formulation of thermal conductivity is similar to Luttinger's 
theory |10l] based on a gravitational field which causes energy or heat current 
to flow. (In his theory, an energy density h{r) is regarded to behave as if it 
had a mass density h(r ) / (? for light velocity c and a very weak gravitational 
field is induced ll|, |l2| by the inhomogeneity of temperature. Of course, 
Luttinger used this idea only conceptually in his formulation.) However, our 
formulation of introducing the above thermal field Et^ and its conjugate heat 
operator Ah through the idea of renormalization of semi-local thermal effects 
is more transparent for deriving the transport coefficient (144 p of the Kubo 
form. Thus we have arrived at a unified "mechanical" formulation of linear 
responses including thermal conductivity ji], lof. 

In particular, to take into account the supply of heat energy is crucial to 
keep the temperature gradient in the thermal conduction. This mechanism 
of energy supply will be formulated in Sections 7 and 8. 

6. Entropy production of thermal disturbance 

We discuss here the entropy production of thermal disturbance using the 
above Kubo formula (jH]). As for the static electric conduction, the entropy 
production for the thermal disturbance is given similarly in the form 

per unit volume and per unit time, as far as the lowest order of is con- 
cerned. In general, we have 

iEr)=<^ (46) 

dt.),J ^' T{t) ^ ' 

for an arbitrary strength of the thermal field. These formulas have been 
derived here on the basis of the above mechanical formulation. 
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In contrast to the electric conduction, the converse argument on the en- 
tropy production, namely to discuss it thermodynamically after the firs- 
principles derivation of it will be very instructive in order to understand 
physically the mechanism of the entropy production for the thermal distur- 
bance. 

The irreversible entropy production per unit volume and per unit time is 
also calculated by 

(f = " FTat) = ^ ^^^^ 

in the lowest order of i?T, where = {AT/L)/T, namely AT = LTE^ for 
the system size L. The heat energy supplied from outside to keep the thermal 
field or thermal gradient (as has been discussed in Section 5) is transformed 
to the relevant system, that is, the energy('Hi)^ = — (Ah ■ -Et)^ is changed 
into the internal energy ('Ho)j, whose increase yields the entropy production 

dS\ 1 d {-Ho), K{E^)El 
dt),^^ T dt T ' 

as in the electric conduction. This gives a unified explanation of the entropy 
production both for the electric conduction and thermal conduction. 

The stationary state for the thermal disturbance can also be formulated [l| 
using the symmetry-separated von Neumann equation, as in the case of the 
electric conduction. 



7. Symmetry-separated von Neumann equation with energy sup- 
ply and new scheme to express stationary states 

In the previous paper [1], we have introduced the symmetry-separated von 
Neumann equation in order to take into account the effect of extracting the 
generated heat outside and to keep the internal energy ("Ho)* constant. In 
the present and succeeding sections, we assume again that the operator "Hi is 
time- independent. However, the energy {T-Li)^ decreases as time t increases. 
In this section, we try to find a new formulation in which both {'Ho)^ and 
{'Hl)^ are time-independent, while there exists the current J = Tijp{t) and 
the entropy production {dS/dt)ij-r is positive, even in steady states. For 
this purpose, we extend here the concept of a steady state by defining the 
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corresponding stationary density- matrix p^^^\t) as 

— = constant operator, (49) 

namely, p^^^^ (t) may contain a linear part with respect to time t as 

p(^^)(t) = p(0)+tp'(0) (50) 

for the density matrix p(0) defined in the limit to — )• — oo. (For t ^ to ^ — oo, 
the above t should be interpreted as te"^ with e — )■ +0 and similarly in 
Eqs.(I71]) and (CSl).) 

The main problem here is to find an extended von Neumann equation to 
describe the process of energy supply without changing the previous physi- 
cal results[i[ on {'Ho)^ (= constant) and {dS/dt)i„{> 0) except (T/i)^ (which 
decreases as time t increases in the previous treatment [l|). Only one addi- 
tional condition of the present new formulation is to require that {'Hi)^ is 
also constant. 

Our strategy to accomplish this scenario is the following: 

(i) First notice that the idea[l[ of symmetry separation of the density ma- 
trix p(t) as p(t) = ps(t) + ps^it) plays a crucial role in formulating energy 
supply from outside to the relevant system as well as in formulating the 
extraction of the generated heat[H. In fact, the effect of energy sup- 
ply is expressed in the partial von Neumann equation concerning the 
antisymmetric part Pa(^) as 

= Uno,pM + 4i^i>Ps(t)] -epa(t) +r/es(t), (51) 
at in in 

where rjesif) denotes the energy source term. 

(ii) The energy source term r]es{t) should be determined so that the solution 
Pa(t) of Eq.ilSTl) coupled with the following equation 

^ = ^[^o,Ps(t)] + ^[^i,Pa(t)] -er(Ps(t) -Po) (52) 

may satisfy the conditions: 
(a) The contribution of the source term ?7es(^) to the derivative d (Hi)^ / dt, 
namely Tr7/i?7es(^) should be non- vanishing. 
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(b) The other quantities {HQ)^ = TiTioPsit), Jp = 1ijps,{t) and 
{dS/dt)i„ should be the same as in the case of no energy supply 
(namely, ?7es(t) = 0). 
Can we find such a desirable operator ?7es(^)? 
(iii) To answer the above question on r]es{t), we try to study what kind of 
the density matrix p{t) satisfies the above two conditions (i) and (ii). 

At first sight, this problem seems difficult. However, we find the following 
favorable relations: 

(a) For any operator function po = Poi'Ho), we have 

TrjiHiPo} = 0; {AB} = ^{AB + BA}, (53) 

and more generally for any operator function p^ = p^('Ho,m) with ?^o,^ = 
"Ho + p'Hi, we have 

TrjWip.} = 0; {A''B} = -^{A''B + A''-^BA + --- + BA'') (54) 

n + 1 

for any positive integer n and for any real value of p, where j = A = 
[A,'Ho]/{ih) and —jF= ['Hi,'Ho]/ {ih). Equation (!53l) corresponds 
to the identity 

Tr ([Hi,Ho]{Hip^}) = Tr ([Hi,HoJ{Hip^}) = 0, (55) 

which can be easily proven using the rule TtABC = TiBCA, and using 
[Hq,^, p^] = 0. Equation can also be proven similarly. 

Note also the following relations 

(b) For p = 0, we have 

Tr(Ho{HiPo}) = 0, (56) 

and more generally 

Tr(Ho{H?"-Vo}) = 0, (57) 
for any positive integer n, as is easily seen from symmetry. 

These relations suggest that the solution of Eq. (l5T|) will take the possible 
form 

Pa(t) = Pa(t) + Pa(t); Pa(t) = (t) {?/lPo} , (58) 
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more generally 

~P.it) = ees(t){?/r^V~o}, (59) 

where Pa.it) is the solution of Eq. llST]) without the source term rj^sit). From 
the condition Pa{to) = 0, we have ees(to) = 0. Note also that the trace of the 
total density matrix 

pit) = pit) + pit) with pit) = pa(t) + Psit) (60) 

is conserved, namely Trp(t) = 1, when Pa,it) takes the form ( 158]) for any 
function of ees(t). Here, ps(t) and /)a(^) are the solutions of the following 
coupled equations without the source term rjesit): 

^Mt) = ^-[•Ho,Ps(t)] + ^[•Hi,pa(t)] - e,(/5, - po), (61) 

^Pa(t) = ^[•Ho,Pa(t)] + ^[Hi,pM - ep.it). (62) 

The tilde density matrices ps and pa denote the contribution from the source 
term rjesit) and they satisfy the following coupled equations: 

= ^[no,Psit)] + ^[^l,Pa(t)], (63) 
^Pa(t) = ^[Ho,Pa(t)] + ^[ni,Psit)]+Vesit). (64) 



Then, what is the source term rjcsit) to give the solution of the form 
(or f l59|) ) in the above coupled equations f l63|) and flMl)? The answer is the 
following: 

Vesit) = ^Mi-Hipo} - ^eUt)['Ho,{'HiPo}] - ^[HupM, (65) 

where 

Ps(t) = ^ e^^°[?/i,pa(s)]e-^««rfs, and pa(t) = eUt){niPo}. 

(66) 

We have only solved inversely Eqs. flBSl) and flM|) . The parameter ees(t) is 
taken to be 

ees(t) = (t + a)e^*e;,(0) (67) 
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for e — )■ +0 and for a constant a satisfying the relation a = ees(0)/egg(0), in 
order to describe a steady state, as will be discussed in the next section. The 
parameter fl67|l has the desired properties that ees{t) is essentially linear in 
time t and ees(— cxd) = 0. 

Thus, we have finally obtained our basic equations to describe our desired 
non-equilibrium phenomena with energy supply from outside to the system 
and with the extraction of generated heat outside. One may say that an 
electric field E in "Hi for electric conduction should be regarded to be an 
external force and consequently that the above energy supply to "Hi may be 
unnecessary. However, it seems reasonable in the case of thermal conductance 
to consider heat energy supply. 



8. True stationary states with both energy supply and heat extrac- 
tion 

First we express our basic equation obtained in Section 7 in the following 
compact form 

^p(t) = -£p(t) + A(t)Po. (68) 
Here, we have introduced the vector p(t) and the super-matrices £ and Csit): 

where the super-operators A, B, C and V are defined jl[ by 

A = l-Ur5no, B = Uj.6ni, C = 1 - uj^o^ = ^e^Un 
<-^r=^, ^e = ^, 6aQ = [A,Q]=AQ-QA. (70) 

Furthermore, the vector Pq and the super-operator ci^'^\t) are defined as 

Po=(|^°) and C^^'\t)po = Vesit). (71) 

(Here we may have put po = Po for simplicity. The unique condition on 
Po is the commutativity of it with T-Lq, namely [Ho,po] = 0.) Then, the 
energy source operator rj^sit) given by Eq. fl65l) with Eq. fl66l) is linear in pq, 
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and consequently the super operator di''^\t) can be defined as above. The 
formal solution of Eq.f l68p is given by 

p(t) = e-(*-*')^£.(t')Porfi' + e-(*-*«)^ {^^^ ■ (72) 

It should be remarked that the above C and C^it) are super matrices whose 
components are super operators such as A, B, C and V. 

In the limit to — oo, Eq. fl72l) can be transformed into the following 
"stationary" solution 

/•oo 

p(«t)(t)= / e-'''C,{t-s)p,ds. (73) 

This can be expressed in the form 

p(^*)(t) = p(0)+tp'(0) (74) 

together with the formula 

A(t)Po = /:s(0)Po + i4(0)Po, (75) 



which are correct for any t as is seen from Eq. (l67I) and Eq.f l71l) with Eq. fl65l) 
for to — oo and e — )■ +0. Here we have 

p(0)= / e-"'C,{-t)p^dt 
Jo 

e'^^CsiO)podt- / e~^^C',{0)pQtdt 
Jo 

= {CMPo - ^"'4(0)Po) , (76) 

and 

poo 

p'(0) = / e-'^'C'MPodt = C-'C'MPo- (77) 
Jo 

Note that Trp'(O) = and consequently that Trp(t) = Trps(t) = 1, as it 
should be. It should be noted that our extended definition of the station- 
ary density matrix fl7^ does not violate the ordinary stationarity of the 
current and Ps{t) except the energy supply proportional to time. (See also 
Appendix for the inverse of the super-matrix C whose elements are 
super-operators.) Equivalently we have 

£p(0) = A(0)Po-p'(0), and Cp'{0) = C'MPo- (78) 
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These relations can also be obtained directly from Eq.f l68p . by expressing 
C^{t) as in Eq.(^E>- 

Thus we have finally arrived at the formal solution of the steady state, 
namely the "stationary" density-matrix p^^^'> = p(0) + tp'(O) with Eqs. fl76|) 
and fl77|) . according to our new definition fH9l) . 

The time-derivative of the energy {^-1)1 is also given by 

j^{'Hi), = TT{'HiVcs{t))-JF-F 

= e'Jt)TT{H,{nM)-JF-F 

= e^,(0)Tr (nlpo) - Jf (79) 

using Eq. fl68|l together with Eqs. fl65l) . fl67j) . fl69|) and fITT]) . and also using the 
relation f l53|) . Here, the nonlinear current Ji;' is given by Eq.f lT3|) . Therefore, 
by imposing the condition 

we arrive at the desired conservation relation of {'Hi){- 

{Hi)^ = constant (81) 

as well as the conservation of ('Ho)^, while the entropy production {dS/dt)-^^^ = 
J F ■ F /T is positive. This yields the true steady state with the irreversible 
(positive) entropy production. 

In more general, the time derivative of the average (Q)^ of any physical 
operator Q = Qs + Qa for the symmetric part and antisymmetric part Qa 
can be calculated by the formula: 

dt dt 
d 



= - (TrQsPs(t) + TrQaPa(t)) 
= jTT{Q.p{t)) 

= -Tr (Q ■ Cp{t)) + Tr (Q ■ C.p^) , (82) 

using the basic equation f l68|) and the vector notation Q defined by = 
(Qs,Qa)- For example, we have 'H^ = {Ho,Hi) and p^{t) = (ps(t), Pa(^)) 
with the inner product ("H • p(t)) = TisPsit) + 'Ha.Pa.it). The above formula 
fl82|) is convenient for studying the average of dynamical derivatives [ij] such 
as the entropy production in steady states, as was shown in Paper I. 
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9. Summary and discussion 



The two types of nonlinear responses which have been formulated here 
look very different from one another, and in fact they give very different ef- 
fects, if we make some approximations in the general expressions (1TB]) and 
These two formulations yield the same physical quantities in different 



forms and consequently we obtain many identities on correlations |25|, |26 



These may be related to fluctuation theorems[27-31], which will be discussed 
somewhere else. In particular, the present Zubarev-type formulation of non- 
linear relaxation will be useful to discuss the relaxation of quantum systems 
to equilibrium[32-36]. 

When we treat the formal solution f ITH]) with Eq. fl2CT]) . we have to be 
careful about the boundedness of 7l(t, to). It is often unbounded below. Even 
the linear approximation of it, that is, Zubarev's expression ( 13T]) has some 
mathematical difficulty. The operator TZi defined by Eq.( l29|) is unbounded, 
because the current operator j is proportional to the momentum operator 
p = Pj in the electric conduction, and because the j is proportional to 
in the thermal conduction, as is seen from Eq. fHTl) . However, T-Lq + TZi 
is bounded below in the electric conduction, while it is still unfounded below 
in the thermal conduction. Thus, the Kubo formula works but Zubarev's 
formula does not in the form (l30l) on the thermal conduction. In any case, 
the perturbational formulas of the transport coefficients up to some finite 
orders of the external force can be used safely, while some compact forms 
including infinite order terms have to be treated carefully in this respect. 

We have formulated the scheme showing the mechanism of energy supply 
to the system from outside like a battery in electric conduction. Thus, we 
have succeeded in constructing a "true steady state" with the constant ('Ho)j 
and ('Hi)j and with positive entropy production, as is shown in FigJTJ The 
symmetrized form {HiPo} in Eq. fISS]) is vital to express energy supply in con- 
trast to the commutator [H, p(t)] in the extended von Neumann equation[37 . 
The present formulation of energy supply will be extended by using which 
is an arbitrary function of the operator "Ho.At = + p'^i introduced above 
Eq.( l54|) . In this extended case, the operator p^ contains both symmetric and 
antisymmetric parts and so does the energy source term r)es{t)- Therefore, the 
situation will become complicated. For example, we have to separate ?7es(^) 
as ries{t) = 'rii'^^\t) + ri^^\t) where rit^\t) (and ri^^\t)) denotes a symmetric 
(and an antisymmetric) part of rjcsit). Correspondingly, the time-evolution 
equation of Ps{t) contains the energy source term ris'^^\t). Then, the energy 
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energy supply 



heat generation 



t 



, dt 



heat extraction 



'dS\ 1 Jf-f 

/irr 



= -Cp{t) + A(t)po 



p^^^>(f)=p(0)+ip'(0) 



Fig. 1: Scheme of irreversible stationary transport phenomena (such as the electric conduc- 
tion) with energy supply and heat extraction. Here, J f denotes the current in a field F, and 
Ps{t) and Pa{t) denote the symmetric and antisymmetric parts of the density matrix p{t), 
respectively. The vector-operators p{t) and £s(^)Po ^'^^ defined by p\t) — {ps{t) , Pait)) 
and (Cs(t)Poy = (crPo, »7cs(^)), respectively. 



(?^o)t changes owing to the effect of r]s (t). This additional energy should 
be extracted outside by increasing the parameter in Eq.( l52i) . in order to 
retain the stationary state as before. Furthermore, we have to be careful 
about the normalization Trp(t) = Trps(^) = 1 in this extended case. Namely 
we have to modify ri!t'^\t) such that Trps(t) = 0. This extended formulation 
describes irreversible transport phenomena sustaining also the energy flow 
without dissipation, through the change of the antisymmetric density-matrix 
Pg.{t). This will give a microscopic or statistical- mechanical explanation of 
"non-equilibrium thermodynamics" [38-42], which will be discussed elsewhere. 

It should also be remarked that the unboundedness of "Hi may be rel- 
evant to the stationary condition (IHOj) . In the case in which the quantity 
Tr ("HiPo) diverges as the system size becomes infinite, we have to make 
infinitesimally small in order to satisfy the condition f l80p . This remark on 
the unboundedness of "Hi (and consequently of Ti = Hq + Tii) is also related 
to the reason why the non-equilibrium operator exp(— /3(?^o + "^i)) with TZi 
defined by Eq. fl29|) works in Zubarev's theory instead of exp(— /3('Ho + ^i)) 
in order to describe transport phenomena. This situation will become more 
clear, if we consider a magnetic system composed of localized spins {Sj}. 
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The Hamiltonian 7i^^^ of this magnetic system given by -HM = -Ho + 
with "Hi = - A^'") ■ H and A^"^^ = /ig '^J is very different from the Hamil- 
tonian of electric conduction 'H'^^^ = + l-Li with T-Li = —A^^^ ■ E and 
A^'^^ = e Vj, in the following three respects: 

(i) T-L^^^ describes an equilibrium state, while Ti^^^ describes a non-equilibrium 
state, (ii) 7{*^™) is bounded, while Ti^^^ is unbounded, and (iii) A*" ^ denotes 



the change of spins, namely spin relaxation 43| (not spin transport), while 



A denotes the electric current, as has been discussed in Paper I. 

One of the remarkable features of the present formulation on stationary 
states is not to use explicit forms of interactions between the system and outer 



systems (outside of Hilbert space[44-46] or two coupled systems |4 71. |48|). but 
to make use of only Tio and Hi with the form Tii = — A - F giving a current 
j = A, as in Kubo's theory on the general frame-work of linear transport 
phenomena. This form of Hi has been shown in the present paper to be valid 
even for thermal conductance. 

It is also very interesting to remark that the time derivative of the entropy 
production <7s(t) = {dSldt)in is derived from Eq. (fT^ as 



j^^s{t) = i jy {j{-zhX)j{t; F)), F\ (83) 



which explains the generalized principle of entropy production in the steady 



state 
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First-principles derivations of dissipation and entropy production for gen- 
eral frequency-dependent external forces such as £^cos(a;t) and H cosiujt) 



will be reported in the near future 49 
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Appendix A. Inverse super-matrix of the super-matrix C whose 
elements are super-operators 

We explain here how to define the inverse super-matrix C^^ of the super- 
matrix of the form 
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where the positions of C and D are different from the ordinary case. Here 
we assume that and exist but that and do not necessarily 
exist. They correspond to A,B,C and V in Eq. flTOj) . namely 



A = e,A, B = -e,B, C = eC, D = -eV. (A.2) 

First we define by the relation CC~^ = 1 (unit super- matrix), namely 
we have 

Y,J^ijiJ^'')jk = 5ik ii,j,k = l,2). (A.3) 
j 

Then, by solving the two-dimensional equation f lA.3p . the inverse super- 
matrix C'^ is easily obtained as 

L, - y^L, - \^_c-iDICA-' C-^DKA-^BC-^ + ' ^ ' 

where /C is defined by /C = (1 — A^^BC^^D)^^ . The operator /C is rewritten 
as /C = (1 — A~^BC^^T>)^^ in the correspondence relation (lA.2p . Finally, 
it should be remarked that the relation C~^C = 1, namely C^^C = LCr^ 
is easily confirmed by explicit calculations with use of be above expressions 
flA.4p or by a mathematical argument. (In fact, any solutions Xx and X2 
of CXx = 1 and X2C = 1 satisfy the relation Xi = 1 x Xi = {X2C)Xi = 
X2{CXi) = X2 X 1 = X2 even for the super- matrices C,Xi and X2 in any 
dimensions.) 
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